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O ■ Abstract 

O 



The paper is devoted to apphcation of recently devised ghost-free Analytic Perturbation Theory 
(APT) for analysis of some QCD observables. We start with the discussion of the main problem of 
^ \ the perturbative QCD — ghost singularities and with the resume of its resolving within the APT. 

w* ' By a few examples in the various energy and momentum transfer regions (with the flavor 

CN . number / = 3, 4 and 5) we demonstrate the effect of improved convergence of the APT modified 

perturbative QCD expansion. 

^ ■ Our first observation is that in the APT analysis the three- loop contribution (~ a^) is as a rule 

CN . numerically inessential. This gives raise a hope for practical solution of the well-known problem of 

z!^ , asymptotic nature of common QFT perturbation series. 

f~^ ' The second result is that a usual perturbative analysis of time-like events with the big vr^ term 

O . in the a^ coefficient is not adequate at s < 2 GeV^ . In particular, this relates to r decay. 

^— ^ \ Then, for the "high" (/ = 5) region it is shown that the common two-loop (NLO, NLLA) pertur- 

■^ ' bation approximation widely used there (at 10 GeV < \/s < 170 GeV) for analysis of shape/events 

Q^' data contains a systematic negative error of a 1-2 per cent level for the extracted as values. 

p ■ Our physical conclusion is that the as{M'^) value averaged over the / = 5 data appreciably 

1^ . differs < as{M'^) >j=5 ~ 0.124 from the currently accepted "world average" (= 0.118). 



'^ : 1 Preamble 



In QCD, a dominant means of theoretical analysis is based on perturbation power expansion 
supported by appropriate renormalization group (RG) summation. This perturbative QCD 
(pQCD) satisfactory correlates the bulk of experimental data in spite of the fact that the RG 
invariant power expansion parameter a^ is not a "small enough" quantity. Nowadays, the 
physically accessible region corresponds to three, four and five (/ = 3, 4, 5) flavor numbers (of 
active quarks). Just in the three-flavor region there lie unphysical singularities of the central 
theoretical object — invariant effective coupling ds ■ These singularities, associated with 
the QCD scale parameter Aj-=3 ~ 400 MeV, complicate theoretical interpretation of data in 
the "small energy" and "small momentum transfer" regions ( y/s, q = \/(f ^ 1 -^ 1-5 GeV). 
On the other hand, as it is well known, their existence contradicts some general statements 
of the local QFT. 

In this paper, we flrst discuss this main problem of the pQCD, the singularities lying in 
the physically accessible domain, and then give a resume of its solution within the recently 
devised ghost-free Analytic Perturbation Theory (APT) that resolves the problem without 
using any additional adjustable parameters. Then, we give some impressive results of the 
APT application for analysis of QCD observables. 
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1.1 Invariant QCD coupling and observables 

Usually, the perturbative QCD part of theoretical contribution to observables in both the 
space- and time-like channels is presented in the form of two- or three-term power expansion 

-— — = l + r(x); r{x) = Cias{x) + C2a1 + C3al + . . . ; x = q^ or = s (1) 

(our coefficients are normalized Ck = Ck t^~^ differently from the commonly adopted Ck , like 
in Refs.|]I|, |^, ^) over powers of the effective QCD coupling a^ which is supposed ad hoc to 
be of the same form as in both the channels, e.g., in the massless three-loop case 



1 
+ 



hi (-ln^L + -ln2L + 21nL-- ) - 86162 lnL+^^ 



(3^L^ [ ^ V 2 2) ' ' 2 

Here, L = ln(x/A^) , and for the beta-function coefficients we use the normalization 
(3{a) = -I3q a^ - /3i a^ - /?2 a^ + . . . = -/?o a^ (l + 61 a + 62 a^ + • • •) , 
that is also free of the vr powers. Numerically, they are of an order of unity 

/?o(/) = ^^~^-^ ; biU) = ^^l~ ^\{^ ; /3o(4 ± 1) = 0.875 ± 0.005 ; 6i(4 ± 1) = 0.490 



(2) 



^^12!^' ^i(/) = 2!r^(33-2/) ' /5o(4±l) = 0.875 ±0.005; 6i(4 ± 1) = 0.490;^?^ 



Meanwhile, the RG notion of invariant coupling was first introduced in QED [Q in the 
space-like region in terms of a real constant z^ of the finite Dyson renormalization trans- 
formation. Just this QED Euclidean invariant charge e(g) is the Fourier transform of the 
space distribution e(r) of electric charge (arising due to vacuum fiuctuations around a point 
"bare" electron) discussed by Dirac [^ in 30s — see Appendix IX in the textbook 0. 

Generally, in the RG formalism (for details, see, e.g., the chapter "Renormalization 
group" in the monograph |^ and/or Section 1 in Ref.|P) the notion of invariant coupling 
g{q) is defined only in the space-like domain. 

In particular, this means that if some observable 0{q^) is physically a function of one 
kinematic Lorentz-invariant space-like argument g^ , then, due to its renormalization in- 
variance, it should be a function of RG invariants only. E.g., in the one-coupling massless 
case 

0(gV/i', 9,) = F {-giqVfi', g,)) with Fig) = 0(1, g) . 

Due to this important property, in the weak coupling case we deal with the functional 
expansion of an observable 0{q^) in powers of invariant coupling g . This is a real foundation 
of QCD power expansion ([l|) in the Euclidean case with x = g^ . At the same time, inside 
the RG formalism, there is no natural means for defining invariant coupling g{s) \ s = —q^ 
and perturbative expansion for an observable 0{s) in the time-like region. 



Nevertheless, in modern practice, people commonly use the same singular expression for 
the QCD effective coupling a^ , like (|^), in both the space- and time-like domains. The 
only price usually paid for this transferring from the Euclidean to Minkowskian region is 
the change of numerical expansion coefficients. The time-like ones Ck>3 = dk — Sk include 
negative "vr^ terms" proportional to vr^ and lower expansion coefficients Ck 



{^Hf)f 



ci , 5i= {-nPof (c2 + -hi ci) . 

o 



(3) 



These (rather essential, as far as n ^^{f = 4 ± 1) = 4.340^;723) structures 5k arise 



12| in the course of analytic continuation from the Euclidean to Minkowskian region. The 



coefficients dk should be treated as genuine /cth-order ones. Just they are calculated via 
the relevant Feynman diagrams. 

Table 1 
Minkowskian Ck and Euclidean rfj = q + Si expansion coefficients and their differences. 



Process 


f 


ci = di 


C2 = d2 


cs 


d3 


^3 


5, 


r decay 


3 


l/vr 


.526 


0.852 


1.389 


0.537 


5.01 


e+e~ 


4 


.318 


.155 


-0.351 


0.111 


0.462 


2.451 


e+e" 


5 


.318 


.143 


-0.413 


-0.023 


0.390 


1.752 


Zq decay 


5 


.318 


.095 


-0.483 


-0.094 


0.390 


1.576 



To demonstrate the importance of the "vr^ terms" , we took the / = 3 case for r-decay, 
the / = 4 , 5 cases for e~^e~ —>■ hadron annihilation and the Zq decay (with / = 5) — 
see Table 1 in which we also give values for the vr^-terms. In the normalization (|1|), all 
coefficients c^ , dk and 6k are of an order of unity. In the / = 4, 5 region the contribution 
63 prevails in C3 and {d^l <^ jcal (see also Table II in Bjorken's review [0). 



1.2 Unphysical singularities 

Let us remind to the reader that the ghost-trouble first discovered in QED in the mid-50s 
(and quite soon in the renormalizable version of pion-nucleon interaction) was considered 
there as a serious argument in favor of inner inconsistency of the whole local QFT. In the 
QED case, the ghost singularity lies far above the mass of the Universe and has no pragmatic 
meaning. 

However, in QCD it lies in the quite physical infrared (IR) region and we are forced to 
face it without any excusing arguments. This means that, if one believes in QCD as in a 
consistent, physically important theory, one has no other possibility as to consider the QCD 
unphysical singularities as an artefact of some approximations used in pQCD. This point of 
view is supported by some lattice simulations and solution of Schwinger-Dyson equations 
e.g.. Section 5.3. in a recent review |]13 . 



see. 



For illustration of the fundamental inconsistency of current pQCD practice connected 
with unphysical singularities, take the well-known relation between the so-called Adler func- 



tion D and the total cross-section ratio i? of a related process 

R{s)ds 



In the case of inclusive e^e^ annihilation into hadrons, R{s) is the ratio of cross-sections 
presented in the form R{s) = 1 + r{s) with a function r expandable in powers of as{s) like 
in Eq.(|l|). At the same time, the Adler function is also used to be presented in the form 
D = 1 + d with d expanded in powers of cxs^q^) ■ 

Here, we face two paradoxes. First, as{q^) and, hence, the perturbative D{q^) obeys 
— see eq.(^ — non-physical singularity at q^ = A^ in evident contradiction with rep- 
resentation (^. Second, the integrand R{s) , being expressed via powers of as{s) , obeys 
non-integrable singularities at s = A^ , which makes the r.h.s. of eq.(^) senseless. 

This second problem is typical of inclusive cross-sections, e.g., for the r hadronic decay. 
Generally, in the current literature it is treated in a very strange way — by shifting the 
contour of integration from the real axis with strong singularities on it into a complex plane. 
However, such a "physical" trick cannot be justified within the theory of complex variable. 

1.3 The "ghost" problem resolving 

Meanwhile, as it is known from the early 80s, the perturbation representation (|^) for the 
Minkowskian observable with the coefficients modified by the vr^-terms is valid only at 
a small parameter vr^/ln (s/A^) values, that is in the region of sufficiently high energies 
W^ = yi > Ae"/2 - 2 GeV . 

Here, it is appropriate to remind the construction devised by Radyushkin ^ and Krasni- 
kov — Pivovarov [|10] (RKP procedure) about 20 years ago. These authors used the integral 
transformation 

Ris) = ^ / -D,,{-z)^R[D{q')] (5) 



2vr _ 
reverse to the Adler relation (^) (that is treated now as integral transformation 



oo 



R{s) ds 



R{s)^D{q-)=q' j^ ^i:^^D{i?(.)} (6) 

for defining modified expansion functions 

Ms) = R[«.'(g')] (7) 

for the perturbative QCD contribution 

r{s) = di2ti(s) + rf22l2(s) + dsQi^is) (8) 

to an observable in the time-like region. 

At the one-loop level, with the effective coupling ai = [Po ln(g^/A^)] one has 

2lS^^ (s) = R \a^}-^] = —— arccos , = — arctan - ; L = In -^ (9) 



and for higher functions 



2l«(.) 



1 



f(i) 



PI [L' + 



TT^ 



%"{s) 



PI [L' + rr'' 



V2-^ 2lJ^(^) 



vrVS 



PI [L' + n^ 



(10) 



which are not powers of 21^ (s) . 

The r.h.s of (pi) at L > can also be presented in the form 



2ti 



(1), 



1 vr 

—— arctan — 
nfjo L 



(9a) 



convenient for the UV analysis. Just this form (9a) was discovered in the early 80s in 
while eqs.(0) in Refs.^] and |]lOl- All these papers dealt with HE behavior 



Refs. M and 

and did not pay proper attention to the region L < . 

On the other hand, expression (P) was first discussed only 15 years later by Milton and 
Solovtsov [15|. Just these authors first made an important observation that expression (P) 
represents a continuous monotone function without unphysical singularity at L = and 
proposed to use it as an effective "Minkowskian QCD coupling" a{s) = 2li(s) in the time- 
like region. 

For the two-loop case, to the popular approximation 



/^oaSop(?^) 



biU) 



In I 



l = ln 



A2 



there corresponds ||, p!6[ 



'^pop\^. 



= 2tf'^°^)(s) 



biL 



L2 + 7r2 



aW(s) - 



bi In [VL^ + 7r2J + 1 



Po 



L2 + 7r2 



(11) 



At L ^ vr, by expanding this expression and 2t2 from (|TOD in powers of vr^/L^ we arrive 
at the 7r2-terms (|^). 

Both the functions (P) and ([Tl|) are monotonically decreasing with a finite IR value 
(5(0) = l//?o(/ = 3) ~ 1.4. They have no singularity at L = 0. Higher functions go to zero 
2lfc(0) = in the IR limit. 

As it has first been noticed in |1^, |18|, by applying the transformation D (^) to functions 
21a;(s), instead of «s(?^) powers, we obtain expressions D[2tfc(s)] = Akiq"^) that are also free 
of unphysical singularities. These functions have been discussed at 90s ||T9[ — ||2^ in the 
context of the so-called "Analytic approach" to perturbative QCD. 

Therefore, this Analytic approach in the Euclidean region and the RKP formulation for 
Minkowskian observables can be united in the single scheme, the "Analytic Perturbation 
Theory" — APT, that has been formulated quite recently in our papers |T^ and |T^. In 
the next Section, we give a short resume of this APT construction and then, in Sections 3 
and 4, present the results of its practical applications. 



2 The APT — a closed theoretical scheme 

The APT scheme closely relates two ghost-free formulations of modified perturbation ex- 
pansion for observables. 

2.1 Relation between Euclidean and Minkowskian 

The first one, that was initiated in the early eighties f^, |T^ and outlined above, changes 



the standard power expansion (|T]) in the time-like region into the nonpower one (|^). It uses 
operation Eq.(|^), that is reverse R = [D]~^ to the one defined by the "Adler relation" @ 
and transforms a real function R{s) of a positive (time-like) argument into a real function 
D{q^) of a positive (space-like) argument. 



By operation R, one can define |T^ the RG-invariant Minkowskian coupling a{s) = 
R [cxs] , and its "effective powers" (0) that are free of ghost singularities. Some examples 
are given by expressions (^, ([T0|) and (|ll|). At the one-loop level, they are related by the 
differential recursion relation A;/5o2t^^ = —{d/dL)Wf^ and are not powers of 2t|^ . 

By applying D to 2tfc(s) , one can "try to return" to the Euclidean domain. However, 
instead of as powers, we arrive at some other functions 

Ak{q^) = D [2tfc] , (12) 

analytic in the cut g^-plane and free of ghost singularities. At the one-loop case 

/?o4^^(^^) = ^r747T^-7^> PI^\^') = v:t7^.+7:P^.. ■■■■ (13) 



These expressions have been originally obtained by other means [|1^, ^ in the mid-90s. 
The first function Ai = aan(Q'^) , an analytic invariant Euclidean coupling, should now be 
treated as a counterpart of the invariant Minkowskian coupling a{s) = 2li(s) . Both a^n and 
a are real monotonically decreasing functions with the same maximum value 

«an(0) = a(0) = l//3o(/ = 3) ~ 1.4 

in the IR liniitQ. 

^ Note that the transition from the usual invariant MS couphng as to the Minkowskian a and Euchdean 
Qfan ones can be understood as a transformation to new renorniahzation schemes. At the one-loop case 

as ^ 5^1) = -- arctan(7r/3oa.) and a, ^ aan^^^ = «« + — f 1 - ei/'^°"») . (14) 

Trpo Po V / 

Here, the first transition looks "quite usual" as a can be expanded in powers of a^ , while the second one in 
the weak coupling case behaves like the identity transformation as far as the second nonperturbative term 
g-i//3oas leaves no "footsteps" in the power expansion. 

For both a^^' and aan the corresponding f3 functions have zero at a = l//3o and are symmetric under 
reflection [a — l/2(3o] — > — [a — l/2/3o] . Moreover, the /3 function for a{s) turns out to be equal to the 
spectral function for aan('Z^) - see below Eq.jl^) at A: = 1 . 



All higher functions vanish Ak{0) = 21^(0) = in this limit. For k > 2 , they oscillate in 
the IR region and form [25, ^ an asymptotic sequence a la Erdelyi. 

The same properties remain valid for a higher-loop case. Explicit expressions for Ak 
and 2lfc at the two-loop case can be written down (see, Refs. [^ and p8|) in terms of a 
special Lambert function. They are illustrated below in Figs la and lb. Note here that to 
relate Euclidean and Minkowskian functions, instead of integral expressions (^ and (^ one 
can use simpler relations, in terms of spectral functions p{a) = QAi^—cr) , 

oo oo 

A(g^/) = -/^^P.(a;/); ^.{s; f) = - [ -pk{a- f) , (15) 

s 

equivalent to expressions Ak{q^) = D [21^] , and 2lfc(s) = R[^fe] • 

Remarkably enough, the mechanism of liberation of unphysical singularities is quite dif- 
ferent. While in the space-like domain it involves nonperturbative, power in g^, structures, in 
the time-like region it is based only upon resummation of the "vr^ terms" . Figuratively, (non- 
perturbative !) analyticization [|TU|, ^, ^ in the g^-channel can be treated as a quantitatively 
distorted reflection (under q"^ —>■ s = — q^) oi (perfectly perturbative) 7r^-resummation in 



the s-channel. This effect of "distorting mirror", first discussed in |jT5[ and p9[, is clearly 
seen in the figures la,b mentioned above. 

This means also that introduction of nonperturbative 1/g^ structures now has got 
another motivation, Eg. 1^1^, independent of the analyticization prescription. 

2.2 Global APT 

In reality, a physical domain includes regions with various "numbers of active quarks", i.e., 
with diverse flavor numbers / = 3, 4, 5 and 6. In each of these regions, we deal with a different 
amount of quark quantum fields, that is with distinct QFT models with corresponding 
Lagrangians. To combine them into a joint picture, the procedure of the threshold matching 
is in use. It establishes relations between renormalization procedures for a model with 
different / values. 

For example, in the MS scheme the matching relation has a simple form 

as{q' = Mj; / - 1) = a.(g' = MJ; /) . (16) 

It defines a "global effective coupling" 

as{q^) = as{q^; f) at M^f_^ < q^ < MJ , 

continuous in the space-like region of positive q^ values with discontinuity of derivatives at 
matching points q^ = M ? . To this global as , there corresponds a discontinuous spectral 
density 

Pfc(a) = p,(a; 3) + J^ ^(a - M]) {p,(a; /) - p^ia; f - 1)} (17) 

/>4 



with pA;(cr; /) = 53q:^'(— o", /) which yields [|T^ |TB| via relations analogous to (jT^) 



Au{q') 



7T J a + q 





da . . „, / X 1 f da . . 



IT J a 

s 



(18) 



the smooth global Euclidean and spline-continuous global Minkowskian expansion functions. 





a aan(0)=i« 


> a (0) 






a'"=350 MeV 


1.0- 


a?Hq';3) I 
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0.5- 


a^iq'^SlJj 


Ka(i)(s;3) 


0.0- 


Q "an 


CP^-^^ s'" (GeV) 

' 1 ' 
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"an/ 














i \ a'"=350 MeV 


0.1 - 
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Figure 1: a - Space-like and time-like global analytic couplings in a few GeV domain with f — 3 
and A^^' = 350 MeV; b - "Distorted mirror symmetry" for global expansion functions. All the curves 
in lb correspond to exact two-loop solutions expressed in terms of Lambert function. 

Here, in Fig. la, by the dotted line we give a usual two-loop effective QCD coupling 
asiq"^) with a singularity at g^ = A^. Meanwhile, the dash-dotted curves represent the one- 
loop APT expressions (^) and ([T3|). The solid APT curves are based on the exact two- loop 
solutions of RG equations and approximate three-loop solutions in the MS scheme. Their 
remarkable coincidence (within the 2-4 per cent) demonstrates reduced sensitivity of the 
APT approach (see, also Refs. [^, |21|, ^) with respect to higher-loop effects in the whole 



Euclidean and Minkowskian regions from IR to UV limits. Fig. lb shows higher two-loop 
functions in comparison with a^n and a powers. 

Generally, functions 21^ and Ak differ from the local ones with a fixed / value. Minkow- 
skian global functions 21^ can be presented via 21a:(s, /) by relations 



«(s) = ais; f) + c(/) ; ^2(3) = ^2(3; f) + C2(/) at 



Mj < s < Mj+i 



(19) 



with shift constants c{f), C2(/) representable via integrals over pk{cr] f -\- n) , n > 1 with 
additional reservations, like c(6) = , related to the asymptotic freedom condition. 

Numerical estimate performed in Ref. |T^ (see also Table 6 in Ref.||2^) for traditional 
values of the QCD scale parameter A3 ~ 300 — 400 MeV 



c(3) ~ 0.02, c(4) ~ 3.10"^ , c(5) ~ 3.10" 



C2(/)^3a(Mj)c(/) 



reveals that these constants are essential in the / = 3, 4 region at a few per cent level for a 
and at ca 10% level for 2I2 . 



Meanwhile, global Eucledean functions Ak{q'^) cannot be related to the local ones Akiq"^, f) 
by simple relations. Nevertheless, numerical calculation shows [27, 28 1 that in the / = 3 
region one has approximately 

«an(g') = aan(g'; 3) + c(3) ; A(g') = A(g'); 3) + C2(3) at Mi<s< Ml . (20) 



3 The APT applications 

3.1 General comments 

In what follows, we abstract ourselves of recent successive use of the Analytic approach to 
hadronic formf actors ||30[| and concentrate on the QCD applications of APT. 

To illustrate a quantitative difference between global APT scheme and common practice 
of data analysis in perturbative QCD, consider a few examples. 

In the usual treatment — see, e.g., Ref. [|I| — the (QCD perturbative part of) Minkow- 
skian observable, like e^e~ annihilation or Zq decay cross-section ratio, is presented in the 
form 



R{s) = i?o (1 + r{s)) ; rpT{s) = Ci a^(s) + C2 a^{s 



C3al{s) + ... 



(21) 

Here, the coefficients Ci , C2 and c^ are not diminishing numerically — see Table 1. A 
rather big negative C3 value comes mainly from the —cin'^f3^/3 term. In the APT, we have 
instead 



TAPTis) = dia{s) + d2 2l2(s) + 4 213(5) + . . . 

Ci 2 and d^ 



(22) 
C3 + Ci7r2/52/3 , the 



with reasonably decreasing Feynman coefficients di^2 
mentioned vr^ term of C3 being "swallowed" by a{s) . 

In the Euclidean channel, instead of power expansion similar to (pT]), we typically have 

dAPT{q^) = dia^q^) + ^2^2(9') + d^A^{q^) +... (23) 

with the same coefficients dk extracted from Feynman diagrams. Here, the modification is 
related to nonperturbative, power in g^, structures like in ([1^). 



Table 2 : Relative contributions (in %) of 1- , 2- and 3-loop terms to observables 



Process 


q or y/s 


f 


PT 


APT 


GLS sum rule 


1.73 GeV 


4 


65 


24 


11 


75 


21 


4 


Bjorken. s.r. 


1.73 GeV 


3 


55 


26 


19 


80 


19 


1 


Inch r-decay 


- 2 GeV 


3 


55 


29 


16 


88 


11 


1 


e^e~ —>■ hadr. 


10 GeV 


4 


96 


8 


-4 


92 


7 


.5 


Zo — > hadr. 


89 GeV 


5 


98.6 


3.7 


-2.3 


96.9 


3.5 


-.4 



In Table 2, we give values of the relative contribution of the first, second, and third terms 
of the r.h.s. in (^TD , (|22|) and (|23|) for the Gross-Llywellin-Smith |^ and Bjorken [^] sum 



9 



rules, T - decay in the vector channel |5^, as well as for e+e" and Zq inclusive cross-sections. 
As it follows from this Table, in the APT case, the three-loop (last) term is very small, and 
being compared with data errors, numerically inessential. This means that, in practice, 

one can use the APT expansions ( P^ and (|23| ) without the last term. 

3.2 Semi— quantitative estimate 

This conclusion can be valuable for the case when the three-loop contribution, i.e., d^ is 
unknown. Here, people use the so-called NLLA approximation, that is common practice in 
the / = 5 region. For the Minkowskian observable, e.g., in the event-shape (see, e.g., Ref. 
|]3^ ) analysis there corresponds the two-term expression 

r(s) =cia,(s) + C2a^(s). (24) 

On the basis of the numerical estimates of Table 1, in such a case, we recommend instead 
to use the two-term APT representation 

t^aIAs) = dia{s) + d2 2t2(s) (25) 

which, at L^ ^ vr^, is equivalent to the three-term expression 

r^(s) = di<as ^'^s f +^2as = Ci a^ + C2a^ - 5sa^, (26) 

i.e., to take into account the known predominant vr^ part of the next coefficient C3 . As it 
follows from the comparison of the last expression with the previous, two-term one (|2^) , 
the ds numerical value extracted from eg. (p6D , for the same measured value Vobs, will differ 
mainly by a positive quantity (e.g., in the / = 5 region with as — 0.12 -^ 0.15) 



[Aa 



7r(5T a 



3 



s 3 



1 + 27id2as 



^-^ 1.225 a^^ 



r-^ 



20-=-100GeV 1 + 0. JU Q^s 



0.002 -i- 0.003 (27) 



that turns out to be numerically important. 

Moreover, in the f = 4: region, where the three-loop (NNLLA) approximation is com- 
monly used in the data analysis, the tt^ term ^4 of the next order turns out also to be 
essential. Hence, we propose there, instead of (|2T|), to use the APT three-term expression 

rfp^{s) = diais) + d2 2l2(s) + ^3 213(5) (28) 

approximately equivalent to the four-term one 

rf'(g) = dias + (J2 al -H C3 a^ - 64 a^ ; C3 = d^ - 63 (29) 

or to 

r^{s) = di<as 5-^as - ^1^ 7r^/3o ^t\ + d2 {a^ - tt'^/Sq a^} + ^3 a^ 

10 



with S3 and 64 defined 0, |I^] in eq.(|^). 

The three- and two-term structures in braces are related to specific expansion functions 
a = 2li and 2I2 defined above (|T3) and entering into the non-power expansion 



To estimate roughly the numerical effect of using this last modified expression (^), we 
take the e"'"e~ inclusive annihilation. For ^/s ^ 3 -^ 5 GeV with a^ ~ 0.28 ^ 0.22 one has 



'Aa 



S)4: 



7164 a^ 
1 + 27r(i2«s 



/=4 



3-=-5GeV 



1.07a 



^,4 



1 + 0.974 a, 



0.005^0.002 



— an important effect on the level of ca 1 -i- 2% . 

Moreover, the {Aas)i correction turns out to be noticeable even in the lower part of the 
/ = 5 region! Indeed, to ^/s ~ 10 -^ 40 GeV with a, ~ 0.20 -7- 0.15 there corresponds 

(Aa,)4|fot40GeV ^ 0.71 at ^ (1-1 - 0.3) ■ 10"^ (< 0. 



3.3 Important warning 

It is essential to note that approximate expressions eqs. (p6D and ( ]29|) are equivalent to the 
exact ones (^) and (^Sj) only in the region L = In (s/A^) ^ vr as it is shown on Fig. 2. 



u.u 

0.5 


\ 

\ 
\ 
\ 

\ 

\ 
\ 


0.4 




0.3 


^ppr /"' -■-■-■-::::::::■ 




Q,s^^'(GeV) 


0.2 


I / 1 1 I I 1 1 1 


1 2 



0.6 



0.8 



L/tt 



Figure 2: Comparison of common QCD coupling cts with the APT global ones ( a , aan ) in the 
q, v^ < 3 GeV region at A3 — 400 MeV . By dash-dotted line we give an approximate Minkowskian 
coupling (BOh. All the curves are taken (see Tables 1,5 and 6 in Ref . p8{ ) for the 3-loop global case. 
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One can see that the curve for approximate Minkowskian couphng 



C^appr i "S 



(s) = as{s) - (7r^/3o73)« 



.^3 



(30) 



that precisely corresponds to the popular approximation (^) (and gives rise to the vr^ term 
in the af coefficient) has a rather peculiar behavior. In the region L > vr it goes rather close 
to the curve for a . For instance, at L ~ tt the relative error of approximation is about 
5 per cent. On the other hand, below L ~ 0.8 vr (i.e., W ~ 1.0 — 1.4 GeV) the distance 
CK — ttappr between curves (error of approximation) increases and at L ~ O.Tvr it blows up 
(better to say "comes down"). 

In particular, at s < 2 GeV^ it is rather desorienting to refer to ds{s) and it is erroneous 
to use aappr{s) and common expansion (^). 

This means that 

below s = 2 GeV^ it is nonadequate to use common as(s) and power expansion eg . (PT|) . 

In other words, we claim that below s = 2 GeV^ it is an intricate business to analyze 
data in terms of the "old good" (but singular) as Q. Here, approximate relation (^) does 
not work as it is illustrated in Fig. 2. 

In this, low-energy Minkowskian/Euclidean region data have to be analyzed in terms of 
nonpower expansion (^2p/(^3p and extracted parameter should be a^n(s) ^a{q^) or A*^^-* . In 
Table 3 we give few numerical examples for the chain 



«« 



,{Mr) ^ &{Mr) ^ A^^^ -^ A(^) ^ as{Mz) 



that allows to study the QCD theoretical compatibility of LE data with the HE ones in the 
APT analysis. 



Table 3 


: Numerica 


chain related LE with HE regions 


a(M.) 


«an(M,) 


A(3) 


A(5) 


as{Mz) 


0.309 


0.332 


450 MeV 


303 MeV 


0.125 


0.292 


0.314 


400 MeV 


260 MeV 


0.121 


0.278 


0.299 


350 MeV 


218 MeV 


0.119 


0.266 


0.286 


300 MeV 


180 MeV 


0.116 



Here, the main element of correlation is the chain A^^^^ ^^ A^^^ ^^ A^^^ that follows from the 
matching condition (|16D 0. 

^In particular, this relates to analysis of r decay. In this connection we would like to attract attention to 
the important paper |33[ that treats the r decay within the APT approach (with effective mass of light quarks 
and the threshold resummation factor) and results in A'^' = 420 MeV that corresponds to aan(M^) = 0.32 
or a{M'^) = 0.30 . At the same time, attempts to interpret results of APT for r decay in terms of as, like, 
e.g., in Ref.|3^, needs some special precaution — see next footnote. A more detailed comment on the r 
decay theoretical analysis will be published elsewhere. 

^Generally, it is possible to use correspondence between aan , 5 and a^ as expressed by relations (p^). 
However, the use of af^{fi'^) at /i < 1 GeV as a QCD parameter could be misleading due to vicinity 
to singularity. For example, at A^^^ = 400 MeV one has as(M^) ~ 0.35 and as{l GeV^) ~ 0.55 to be 
compared with a^niM^) ^ 0.31 and aan(l GeV^) ~ 0.40. 
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4 Quantitative illustration 

Consider now a few cases in the / = 5 region. 

T decay. According to the Particle Data Group (PDG) overview (see their Fig. 9.1 on 
page 88 of Ref.!!]]), this is (with as(M|) ~ 0.170 and as(M|) = 0.114) one of the most 
"annoying" points of their summary of as(M|) values. It is also singled out theoretically. 
The expression for the ratio of decay widths starts with the cubic termQ 

R{T)=Roal{^M^){l-eias) with ^ < 0.5 and Ci(0 ^ 1 • 

Due to this, the vr^ corrections^ corresponding to the APT expression 

/?apt(T) = i?o2l3(eM2)(l - ei2t4) (31) 

are rather big Sts ~ af (1 — 2(7r/3o)^Q;^) , 2I4 ~ a^ [1 — {10/ 3) {it (3o)'^al] in the region with 
7r^/?o(5) = 3.57 and as{^My) ~ 0.2 . As a crude estimate (taken from a^ — > 2I3 only), 

Aas{M^) = ^ (vr/?o)' «'(M|) ~ 0.0123 , 

which corresponds to 

Aa,(M|) = 0.006 with resulting a,(M|) = 0.120 . (32) 

One should note here, that this estimate is rather crude and gives only indication of the 
order of magnitude. 

The NNLO case. Now, let us turn to a few cases analyzed by the three-term expan- 
sion formula ([^). For the first example, take e~^e~ hadron annihilation at y/s = 42 GeV and 
llGeV. 

A common form (see, e.g., Eq.(15) in Ref. 0) of theoretical presentation of the QCD 
correction in our normalization looks like 

r,+e- ( v^) = 0.318a,(s) + 0.143 a^ - 0.413 a^ . 

In the standard PT analysis, one has (see, e.g.. Table 3) q;s(42^) = 0.144 that corresponds 
to re+e-(42) ^ 0.0476 . Along with the APT prescription, one should use 

Te+e- iVs) = 0.318 &{s) + 0.143 2l2(s) - 0.023 ^^{s) , (33) 

which yields a{A2^) = 0.142 -^ a,(422) = 0.145 and a,(M|) = 0.127 to be compared with 
q;s(M|) = 0.126 under a usual analysis. 

Quite analogously, with as(ll^) = 0.200 and re+e-(ll) — 0.0661 we obtain via (^) 
a(ll2) = 0.190 that corresponds to as{M^) = 0.129 instead of 0.130. 

''See, e.g., eq.(9.16) in Ref.§. 

^First proposal of taking into account this effect in the T decay was discussed [Q more than a quarter 
of century ago. Nevertheless, in current practice it is completely forgotten. 
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For the next example, take the Z^ inclusive decay. The observed ratio 



Rz = r(Zo -^ hadrons)/T{Zo -> leptons) = 20.783 ± .029 

can be written down as follows: Rz = -Rq (1 + '"z(^l)) with Rq = 19.93 . A common form 
(see, e.g., Eq.(15) in Ref.p|) of presenting the QCD correction rz looks like 



rziMz 



-2 



0.3326a, + 0.0952 at - 0.483 a 



^,3 



To [?"z]o;,s ~ 0.04184 there corresponds as{Mz) 
APT case, from 



0.124 with A. 



(5^ 
MS 



292 MeV . In the 



APT 



{Mz) = 0.3326 a(M|) +0.0952 2t2(M|) -0.094 2l3(M|) 



(34) 



we obtain a(M|) = 0.122 and q;s(M|) = 0.124. Note that here the three-term approxima- 
tion (Ij) gives the same relation between the cXsiM"^) and a(M|) values. 

Nevertheless, in accordance with our preliminary estimate for the {/\as)^ role, even the 
so-called NNLO theory needs some vr^ correction in the W = y/s < 50 GeV region. 

The NLO case. Now, turn to those experiments in the HE (/ = 5) Minkowskian region 
(mainly with a shape analysis) that usually are confronted with the two-term expression 
(|3). As it has been shown above (^Tf), the main theoretical error here can be expressed in 



Table 4 : The APT revised|] part (/ = 5) of Bethke's Table 6 





v^ 


loops 


tts (s) 


asifnl) 


a, (s) 


as{ml) 


Process 


GeV 


No 


ref.[2] 


ref.[2] 


APT 


APT 


T- decay ] 


9.5 


2 


.170 


.114 


.182 


.120 (+6) 


e+e"[cr/jarf] 


10.5 


3 


.200 


.130 


.198 


.129(-1) 


e~^e~[j & sh] 


22.0 


2 


.161 


.124 


.166 


.127(+3) 


e~^e~[j L sh] 


35.0 


2 


.145 


.123 


.149 


.126(+3) 


e^e-[ahad] 


42.4 


3 


.144 


.126 


.145 


.127(+1) 


e~^e~[j Szsh] 


44.0 


2 


.139 


.123 


.142 


.126(+3) 


e~^e~[j &:sh] 


58 


2 


.132 


.123 


.135 


.125(+2) 


Zq -^ had. 


91.2 


3 


.124 


.124 


.124 


.124(0) 


e~^e~[j L sh] 


91.2 


2 


.121 


.121 


.123 


.123(+2) 


V 




2 








...(+2) 


e+e~[j &s/i] 


189 


2 


.110 


.123 


.112 


.125(+2) 



Averaged < as{M'^) >/=5 values 



0.121; 



0.124 



""j & sh" — jets and shapes; Figures in brackets in the last column give the 
difference Aq:s(M^) between common and APT values. 
''Taken from Ref.||. 

the form 



fAa, 



|/=5 
l20-M00GeV 



1.225 a 



;^,3/ 



0.002 + 0.003. 



(35) 
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An adequate expression for the equivalent shift of the as(M|) value is 

[Aa,(M|)]3 = 1.225a,(s)«,(M|)2. 



(36) 



We give the results of our approximate APT calculations, mainly by Eqs. (|35|) and (|36|), 
in the form of Table 4 and Figure 3. In the last column of Table 4 in brackets, we indicate 
the difference between the APT and usual analysis. The results of the three-loop analysis 
are marked by bold figures. Dots in the lower part of the Table correspond to shape-events 
data for energies W = 133, 161, 172 and 183 GeV with the same positive shift 0.002 for the 
the extracted a, values. 



a^.(s) 



0.18 



0.16 _ 



0.14 _ 



0.12 



270 0.1235 

213 0.1184 




Figure 3: The APT analysis for 

by the tt^ correction ( p5[ 

dot-and-dash curve, we give the standard a 

from Fig. 10 of paper H. 



Sohd APT curve relates to AJ^^g 



5 time-like region. Crosses (-f ) differ from circles 

270MeVand as(Af|) == 0.124 . By 

(at A(5) ='213MeV and as(M|) = 0.118) taken 



In Fig. 3, by open and hatched circles we give two- and three-loops data from Fig. 10 of 
paper 0. The only exclusion is the T decay taken from Table X of the same paper. By 
crosses, we marked the "APT values" calculated approximately by Eq. (|35D . 

For clearness of the vr^ effect, we skipped the error bars. They are the same as in the 
mentioned Bethke's figure and we used them for calculating x^ ■ 

Let us note that our average 0.121 over events from Table 6 of Bethke's review nicely 
correlates with recent data of the same author (see Summary of Ref.[3G]). The best x^ fit 
yields a.(M|)[2] = 0.1214 andQ 



as{Ml)APT 



0.1235, 



^This value, corresponding to A*^^) — 290 MeV, is supported by recent analysis ||33 
A(3) = 420 MeV ; compare with Table 3. 



of T decay that gives 
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This new Xapt i^ smaller Xj^ptI'Xj^t — O-''^^ than the usual one. This illustrates the effec- 
tiveness of the APT procedure in the region far enough from the ghost singularity. 

5 Conclusion 

It is a common standpoint that in QCD it is legitimate to use the power in a^ expansion for 
observables in the low energy (low momentum transfer) region. At the same time, there exist 
rather general (and old |^) arguments in favor of nonanalyticity of the S matrix elements 



at the origin |5^ of the complex plane of the expansion parameter a variable. This, in turn, 
implies that common perturbation expansion has no domain of convergence. Technically, 
this corresponds to the factorial growth (~ n! ) of expansion coefficients (like dn or Vn ) at 
large n |3^, ^. In QCD, with its "not small enough" as values in the region below 10 GeV 
it is a popular belief that one does face an asymptotic nature of perturbation expansion by 
observing approximate equality of relative contributions of the second (a^) and the third 
(Og) terms into observable, like in all PT columns of Table 2. 

Our first qualitative result consists in observation that convergence properties of the APT 
expansions drastically differ from the usual PT ones. 

The evidently better practical convergence of the APT series for the Euclidean observable, 
as it has been demonstrated in the right part of Table 2, probably means that essential 
singularity at Os = is adequately taken into account by new expansion functions Ak{q^) ■ 
On the other hand, in the time-like region the improved approximation property of the 
APT expansion over 2tfc(s) has a bit different nature, being related, in our opinion, to the 
non-uniform convergence of the standard PT expansion for Minkowskian observables. In 
any case, from a practical point of view: 

1. In the APT, one can use the nonpower expansions ( P^ and ( ^3]) without the last term. 



The next point, discussed in Section |3^, refers to a more specific issue connected with 
current practice of the Minkowskian observable analysis in the low-energy (s < 3 GeV^) 
region (like, e.g., inclusive r decay). As it has been shown - see Fig. 2 — 

2. Below 2 GeV^ it is impossible to use the common power expansion (|l|) for a time-like 
observable. 

Second group of our results is of a quantitative nature: 

3. Effective positive shift Aa^ ^ +0.002 in the upper half (> 50 GeV) of the / = 5 region 
for all time-like events that have been analyzed up to now in the NLO mode. 

4. Effective shift Acts > +0.003 in the lower half (10 + 50 GeV) of the / = 5 region for all 
time-like events that have been analyzed in the NLO mode. 

5. The new value 

a,(M|) = 0.124, (37) 

obtained by averaging new APT results over the / = 5 region. 
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The quantitative results are based on the new APT nonpower expansion (^ and plausible 
hypothesis on the 7r^-term prevalence in common expansion coefficients for observables in 
the Minkowskian domain. The hypothesis has some preliminary support — see Table 1 — 
but needs to be checked in more detail. 

Nevertheless, our result (|37D being taken as granted raises two physical questions: 

- The issue of self-consistency of QCD invariant coupling behavior between the "medium 
(/ = 3,4)" and "high (/ = 5,6)" regions. 

Here, the more detailed APT analysis of data on DIS, heavy quarkonium decays and some 
other processes are in order. As it has been mentioned above, fresh APT analysis of r-decay 



33| seems to support such a correlation with A3 ~ 400 -=- 450 MeV and A5 ~ 290 MeV. 

- The new "enlarged value" (^) can influence various physical speculations in the very 
HE region, in particular on superpartner masses in MSSM GUT constructions - compare, 
e.g., with recent attempts |JT| in this direction. 
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